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Transportation of ions and neutrons in matter is of direct interest in several 
technologically important and scientific areas, including space radiation, cosmic ray 
propagation studies in galactic medium, nuclear power plants and radiological effects that 
impact industrial and public health. For the proper assessment of radiation exposure, both 
reliable transport codes and accurate data are needed. Nuclear cross section data is one of 
the essential inputs into the transport codes. In order to obtain an accurate parametrization 
of cross section data, theoretical input is indispensable especially for processes where there 
is little or no experimental data available. 

In this grant period work has been done on the studies of the use of relativistic 
equations and their one-body limits. The results will be useful in choosing appropriate 
effective one-body equation for reaction calculations. Work has also been done to improve 
upon the data base needed for the transport codes used in the studies of radiation transport 
and shielding for space exploration and high speed flight transportation. A 
phenomenological model was developed for the total absoiption cross sections valid for 
any system of charged and/or uncharged collision pairs for the entire energy range. The 
success of the model is gratifying. It is being used by other federal agencies, national labs 
and universities. A list of publications based on the work during the grant period is given 
below and copies are enclosed with this report. 
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ONE-BODY LIMITS OF TWO-BODY RELATIVISTIC 
EQUATIONS 


Trina L. Christian and Khin Maung Maung 
Department of Physics , Hampton University, Hampton VA 23668 

John W. Norbury 

Department of Physics, University of Wisconsin , Milwaukee , WI 54601 

The bound state spectra of the Blankenbecler- Sugar, Kadyshevsky and Gross 
equations with a linearly confining interaction are studied. As a result of the lack 
of retardation in the interaction, all of these equations approach the same one- 
body limit, but the rates of approach are different. Numerical results show that 
the one-body limit is reached when 7712/m i = 10 for this particular interaction. 


In a two particle system of masses mi and m 2 , in the limit that mi or m 2 
goes to infinity the relativistic equation for this system should be one which 
describes the motion of the lighter particle in the instantaneous field of the 
heavier particle; the relative energy between the particles is zero. This is known 
as the one-body limit. The Blankenbecler-Sugar 1 , Kadyshevsky 2 , and Gross 3 
equations are studied to determine how the one-body limit is approached as m 2 
is increased with mi fixed. The momentum space two-body relativistic bound 
state equations are solved with a confining potential without retardation. 4 

A wide class of relativistic two-body propagators can be written as 

g(s',s) = -2m [ -^-^—6 {+) (A)6 {+) (B)ds' 

J s' - s + irj 

where k and P are the relative and total four-momenta and /(s',s) is any 
function that satisfies the constraint f(s,s) = 1. Therefore in principle there 
exist an infinite number of three-dimensional equations. 5 The choices for A and 
B in the arguments of the delta functions are chosen to distinguish between 
one particle on-mass-shell systems and both particles equally off- mass-shell. 

'' For one particle on-mass-shell systems, Type I, A = (aP + k) 2 — m 2 and 
B = (P' -aP- k) 2 - m\ and A = ( aP ' + k) - m 2 and B = ( bP ' - k ) 2 - m 2 
for Type II, both particles equally off- mass-shell. In order to have a one-body 
limit in three-dimensional equations, we use the Wightman-Garding variable^ • 
7 , a and 6; a = (s + m 2 - m 2 )/2s and 6 = (s + m\ — m 2 )/2s where s is the total 
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four-momentum squared. Explicit forms of propagators in the center-of-mass 
frame have been tabulated 8 The Blankenbecler-Sugar propagator is Type II, 
while the Kadyshevsky and Gross propagators are Type I. 

In these preliminary calculations, we study the bound state spectrum of 
these equations with a linear confining interaction without retardation. The 
figures represent each of the three equations evaluated with increasing m 2 . In 
figures (1) to (4) the mass ratios are 1,2,5 and 10 respectively. In each figure, 
squares, diamonds and triangles represent the Kadyshevsky, Blankenbecler- 
Sugar and Gross results. 


m 2/m 1 ■ 1 



m 2/m 1 a 2 
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The curves display important revelations: for this particular interaction, 
the various equations approach the same one-body limit at different rates, and 
a mass ratio is determined for a two particle bound state system to reach the 
one-body limit at m 2 to mi equal 10. 

The results can be understood in the following way. As mi or m 2 goes to 
infinity the choice of type becomes insignificant. The function /(s',s) goes to 
unity as well if one particle mass goes to infinity. Therefore for any nonretarded 
interaction, they all have the same one-body limit. Although a one-body limit 
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was reached, the Klein-Gordon result was not obtained. This shortcoming is 
attributed to additional multiplicative terms that vary from one propagator to 
the next but nevertheless prevent the Klein-Gordon limit from being reached. 1 2 3 4 5 6 7 8 9 
However, /(s', s) does not equal one in general. This accounts for the difference 
in the rates of approach to the one-body limit for the different equations. 

In conclusion, it is shown that the same one-body limit is reached for Type 
I and Type II equations without retardation in the interaction. This can be 
helpful when evaluating a two particle system; the results suggest that a one- 
body equation should not be used if the mass ratios of the particles are not 
consistent with the results obtained above, at least for this particular inter- 
action. This type of study should also be considered for scattering problems 
where the use of a one-body equation is very common. 
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Abstract 

We present a simple universal parameterization of total reaction cross sections for any system of colliding nuclei valid for 
the entire energy range from a few A MeV to a few A GeV. The universal picture presented here treats the proton-nucleus 
collision as a special case of the nucleus-nucleus collision, where the projectile has charge and mass number one. The 
parameters are associated with the physics of the collision system. In general terms Coulomb interaction modifies cross 
sections at lower energies and the effects of Pauli blocking are important at higher energies. The agreement between the 
calculated and experimental data is better than all earlier published results. 


1. Introduction 

The transportation of energetic ions in bulk matter is 
of direct interest in several areas [ 1 ] including shielding 
against ions originating from either space radiations or ter- 
restrial accelerators, cosmic ray propagation studies in galac- 
tic medium or radiobiological effects resulting from the work 
place or clinical exposures. For carcinogenesis, terrestrial ra- 
diation therapy, and radiobiological research, knowledge of 
the beam composition and interactions is necessary to prop- 
erly evaluate the effects on human and animal tissues. For 
the proper assessment of radiation exposures both reliable 
transport codes and accurate input parameters are needed. 

One such important input is the total reaction cross sec- 
tion, defined as the total minus the elastic cross sections for 
two colliding ions: 

CTR = CTJ — (T el . (1) 

In view of its importance the total reaction cross section has 
been extensively studied both theoretically [1-14] and ex- 
perimentally [ 15-20] for the past five decades. A detailed 
list of references is given in Refs. [1,13,16]. Empirical pre- 
scriptions have been developed [2-4,10,1 1,13] for the total 
reaction cross sections working in various energy ranges and 
combination of interacting ions. The present model works in 
all energy ranges for any combination of interacting ions in- 
cluding proton-nucleus collisions and is more accurate than 
earlier reported empirical models. 


■“Corresponding author. Tel. +1 804 864 1467, fax +1 804 864 8094, 
e-mail rkt@hesbl.larc.nasa.gov 


2. Model description 

The present model is an updated and revised version of 
the empirical model developed at NASA Langley Research 
Center and reported earlier [10]. Most of the empirical mod- 
els approximate the total reaction cross section of the Bradt- 
Peters form: 

crabs = 7rro 2 (Ap /3 + Ay 3 - <5) 2 , (2) 

where ro is energy independent and 8 is either an energy- 
independent or energy dependent parameter, and Ap and At 
are the projectile and target mass numbers, respectively. This 
form of parameterization works nicely for higher energies. 
However, for lower energies Coulomb interaction becomes 
important and modifies the reaction cross sections signif- 
icantly. In addition, strong absorption models suggest an 
energy dependence of the interaction radius. Incorporating 
these effects, and other effects discussed later in the text, we 
propose the following form for the reaction cross section: 

cr R = -n rl(Af + Aj /3 + Se) 2 { 1 - 5/Ecm). (3) 

We notice that the Coulomb interaction, where ro = 
1.1 fm, and E cm is in MeV, modifies the cross sections at 
lower energies and gets less important as the energy in- 
creases (typically after several tens of A MeV). In Eq. (3) 
B is the energy dependent Coulomb interaction barrier 
(right hand factor in Eq. (3)), and is given by 

B= 1.44Z p Z t //?, (4) 

where 

R = rp + ry + 1 .2(Ap 3 + a| /3 )/£^ 3 , (5) 


0168-583X/96/$15.00 Copyright (c) 1996 Elsevier Science B.V. All rights reserved 
/>// SO 1 68-583 X( 96) 0033 1-X 


348 


R.K. Tripat hi et al./Nucl. Instr. and Meth. in Phys. Res. B 1 17 (1996) 347-349 


P + 9 4 Be 



Fig. 1 . Absorption cross sections for the proton-beryllium collision as a 
function of proton kinetic energy. The solid line represents present model 
and the experimental data are from Ref. ( 1 5 1 . 

with (i = P, T) 

17 = 1.29(r/)rms. (6) 

There is an energy dependence in the reaction cross sec- 
tion at intermediate and higher energies mainly due to two 
effects - transparency and Pauli blocking. This is taken into 
account in 8e , which is given by 

S E = 1 .855 + 0. 1 6-4 t ~ C E + 0.91 ( — ~ 2Zt) Zp , (7) 

E'L 3 *tA, 

where S is the mass asymmetry term and is given by 

S = a'^a'^Ka'J 3 + a}/ 3 ) (8) 

and is related to the volume overlap of the collision system. 
The last term on the right hand side of Eq. (7) accounts for 
the isotope dependence of the reaction cross section. The 
term Ce is related to the transparency and Pauli blocking 
and is given by 

C e = D( 1 -exp(— £/40)) - 0.292 exp(-£/792) 

xcos( 0.229 fi"' 453 ). (9) 

Here D is related to the density dependence of the colliding 
system scaled with respect to the density of the C+C system, 
i.e.: 

D = 1 J5(pA P + PAi)/{pA c + pA C ) • (10) 

The density of a nucleus is calculated in the hard sphere 
model [24], and for a nucleus of mass number A-, is given 
by 

pA, = Ail J-7T r 3 , (11) 

where the radius of the nucleus r, is defined in Eq. (6) 
with the root-mean-square radius, (/*,•) rms, obtained directly 
from experiment [25]. There is interesting physics associ- 
ated with the constant D. This in effect simulates the modi- 
fications of the reaction cross sections due to Pauli blocking. 


P + 27 i 3 AI 



Fig. 2. Absorption cross sections for the proton-aluminum collision as a 
function of proton kinetic energy. The solid line represents the present 
model and the experimental data are from Ref. [15]. 

This effect is new and has not been taken into account in 
other empirical calculations. This helps present a universal 
picture of the reaction cross sections. 

At lower energies (below several tens of A MeV) where 
the overlap of interacting nuclei is small (and where 
Coulomb interaction modifies the reaction cross sections 
significantly) the modifications of the cross sections due to 
Pauli blocking are small, and gradually play an increasing 
role as the energy increases, since this leads to higher den- 
sities where Pauli blocking gets increasingly important. In- 
terestingly enough for the proton-nucleus case, since there 
is not much compression effect, a single constant value of 
D = 2.05 gives very good results for all proton-nucleus 
collisions. For alpha-nucleus collisions, where there is a 
little compression, the best value of D is given by 

D = 2.77 - 8.0 x 10 -3 At+ 1.8 x \0~ 5 a\ 

— 0.8/( 1 +exp(250— £)/75). (12) 

For lithium nuclei because of the “halos”, compression is 
less and hence the Pauli blocking effect is less important and 
a reduced value of D/3 gives better results for the reaction 
cross sections at the intermediate and higher energies. 

There are no adjustable parameters in the model except 
that for proton-nucleus collisions this method of calculat- 
ing the Coulomb energy underestimates its value for the 
very light closed shell nuclei of alpha and carbon, and these 
should be increased by a factor of 27 and 3.5 respectively 
for a better fit. 

3. Results/conclusions 

Typical results obtained from the model are shown in 
Figs. 1 through 5. Agreement with experimental data is ex- 
cellent and is better than all other empirical models reported 
earlier. This is particularly important in view of the fact 
that the agreement is excellent throughout the whole energy 
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Fig. 3. Absorption cross sections for the alpha-carbon collision as a function 
of incident ion kinetic energy. The solid line represents present model and 
the experimental data are from Refs. | 21,22). 
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Fig. 4. Absorption cross sections for the carbon-carbon collision as a 
function of incident ion kinetic energy. The solid line represents the present 
model and the experimental data are from Refs. [ 1 ,9, 13.16]. 
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Fig. 5. Absorption cross sections for the calcium-calcium collision as a 
function of incident ion kinetic energy. The solid line represents the present 
model and the experimental data are from Refs. [1,9,13,16]. 


range - up to a few A GeV. The model has been tested with 
all the available data for projectiles proton through krypton 
and targets alpha through bismuth for the energy range from 
a few A MeV up to a few A GeV and is found to give excel- 
lent results for all the systems throughout the energy range. 
In view of the simplicity and accuracy of the model it is a 
welcome improvement for transport calculations. 

It will be interesting to see how the model compares with 
the new experimental data as and when these become avail- 
able. 
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Abstract 

A recent parameterization (here after referred as paper I, Ref. [4]) of absorption cross sections for any system of charged 
ions collisions including proton -nucleus collisions, is extended for neutron-nucleus collisions valid from ~ 1 MeV to a 
few GeV, thus providing a comprehensive picture of absorption cross sections for any system of collision pair (charged 
and/or uncharged). The parameters are associated with the physics of the problem. At lower energies, the optical potential 
at the surface is important and the Pauli operator plays an increasingly important role at intermediate energies. The 
agreement between the calculated and experimental data is better than earlier published results. 


1. Introduction 

The transport of neutrons in matter is of direct interest 
in several technologically important and scientific areas 
[1-3] including space radiations, cosmic ray propagation 
studies in galactic medium, nuclear power plants, radiobio- 
logical effects impacting on industrial and public health. 
For the proper assessment to radiation exposures both 
reliable transport codes and accurate input data are needed. 

An important ingredient of the input data is the total 
absorption (reaction) cross section, defined as the total 
minus the elastic cross sections for two colliding ions: 

<r R = <T r - (T et . ( 1 ) 

Recently, we have developed a simple accurate formalism 
(paper I and Ref. [5]) for the total absorption cross sections 
for any system of charged colliding ions including protons. 
The present work extends the formalism to neutron-nucleus 
collisions, thus presenting a simple accurate comprehen- 
sive formalism for the absorption cross sections for any 
combination of colliding ions (charged and/or uncharged) 
valid for the entire energy range. We have retained the 
high energy features of the model [6] developed here at 
NASA Langley Research Center (LaRC) where it is accu- 
rate, and improved on its behavior at lower energies within 
the universal formalism recently developed (paper I and 
Ref. [5]). The other commonly used model [7] gives good 
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results for higher energies but shows increasingly larger 
deviations from experiment for lower energies. 


2. Model description 

Most of the empirical models approximate total reac- 
tion cross section for charged ion collisions of Bradt-Peters 
form [8], but have very different form for the neutron- 
nucleus collisions. We have maintained the uniform con- 
sistent picture for all the systems proton-nucleus, nu- 
cleus-nucleus, and now for the neutron-nucleus colli- 
sions, and have written total reaction cross sections for 
neutron-nucleus collisions also of the Bradt-Peters form: 

Or” nr^A'/ 3 +A 1 -/* - S)\ (2) 

where r 0 is energy independent and 8 is either energy-in- 
dependent or dependent parameter, and A ? and A T are the 
projectile and target mass numbers respectively. 

It is well known that this type of parameterization is 
suited for higher energies and is modified through Coulomb 
interaction at lower energies (paper I and Ref. [5]) for the 
charged ion collisions. For the neutron-nucleus collisions 
there is no Coulomb interaction, but the total reaction cross 
section, in this case, is modified by the strength of the 
imaginary part of the optical potential at the surface. Since 
we are using this form of parameterization for the neu- 
tron-nucleus case for the first time (which helps to pro- 
vide the unified consistent and accurate picture for the 
total reaction cross sections for any system of colliding 
nuclei for the entire energy range), we introduce a low 
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energy multiplier (X m ) accounting for the strength of the 
optical model interaction. In addition, the effects of the 
transparency and Pauli blocking are taken into account 
through the energy dependent parameter, S E , defined latter 
in the text. We, therefore, propose the following similar 
form for the reaction cross sections for the neutron- 
nucleus: 


The density of a nucleus is calculated in the hard sphere 
model, and for a nucleus A x is given by 

47 r 

Pa,~ *,/—*?, (II) 

where r x is the equivalent sphere radius and is related to 
the r mis j radius by 


1 !rr o( /l P / 3 + ' 4 T /3 + 5 E)‘^m. ( 3 ) 

where r 0 = 1.1 fm. The low energy multiplier X m is given 
by 

(i) For A t < 200 

<4 > 

where S L = 0.6 for A T < 12, S L = 1.6 for A T = 12, S L = 1 
otherwise with 


X, = 2.83 — 3.1 X 10 " 2 Aj + 1 .7 X 10 ~W T . (5) 

For the low energy, shell effects play an important role in 
light nuclei and the optical potential at the surface varies 
differently with these effects showing variation in S L . 

(ii) For A r > 200 


X 


rn 


- 0.3exp 



[ I — exp(0.9 — (6) 


The effects of the energy dependence at intermediate 
and higher energies due to Pauli blocking and transparency 
are taken into account in 5 E , which is defined as, 

S E = 1.855+0. 16 S/E' c Q 

— C E + 0.9 1 ( A r — 2 Z t )Z p /(A t A p ), (7) 

where S is the mass asymmetry term and is given by 


a'/ 3 aY } 
A'j' + Ay* 


( 8 ) 


and is related to the volume overlap of the collision 
system. For the results reported here A ? = 1. The last term 
on the right hand side of Eq. (7) accounts for the isotope 
dependence of the reaction cross section. The term C E is 
related to the transparency and Pauli blocking and is given 
by 

C e = D[| - exp( -£/r,)] - 0.292 exp(-£/792) 


X cos(0.229£°' 453 ), 


( 9 ) 


where 7’, = 40 (except for 11 < A T < 40 a value of 7', = 30 
is recommended for a better fit). Here D is related to the 
density dependence of the colliding system and the global 
value is given by 


0.538 
Pa v + Pa, 


( 10 ) 


'•i=l-29r nra>i , (12) 

where r mis ; values are taken from the experiment [9], 

There is interesting physics associated with the constant 
D. This in effect simulates the modifications of the reac- 
tion cross sections due to Pauli blocking. This effect is 
new and was introduced in our charged particles reaction 
cross section work (paper I), and is being used for the first 
time in neutron-nucleus work. This is an important physi- 
cal ingredient of the model and helps present a unified 
picture of the reaction cross sections for any system of 
colliding particles. 

Let us elaborate on the choice of D in Eq. (9). This has 
contributions from the global characteristics of nuclei and 
also takes into account the specific characteristics related 
to their stability. The density of a nucleus is a very useful 
global property. We, therefore, associate the global value 
D g to the nuclear densities (Eq. (11)) of the colliding 
system. Consequently, for most systems, we have: 

D = D g , (13) 

where D g has been defined in Eq. (10). This value of D is 
modified in certain groups of nuclei to account for their 
special features in the following way: 

The shell structure plays an important role in the size of 
a light nucleus and hence its cross sections. These effects 
are more prevalent upto double closed sd shell nuclei 
U<40) and manifest themselves as a tightly bound 
nucleus and also as isotopic effects. Hence for this group 
of nuclei, by taking into account these special character- 
istics, the best value of D is given by: 

For A T < 40 

£> = D g - 1 .5( A T — 2Z t )/A t 

+ 0.25/{ 1 4- exp[( E — 170)/10()]}. (14) 

The second term on the RHS of Eq. (14) accounts for 
the isotopic effects and the last energy dependent term 
accounts for the shell effects upto doubly closed sd shell 
nuclei. For still heavier nuclei shell effects start becoming 
less important but the isotopic effects persist on. We, 
therefore, found better agreement using the following ex- 
pression: 

For 40 < A T < 60 

D = Dg— 1 .5( A T — 2Z t )/A t . (15) 

For even higher nuclei (/l > 60) special characteristics 
do not dominate the stability of a nucleus and the global 
value of D defined in Eq. (13) works well. However, for 
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Fig. 1. Total absorption cross sections for n + 4 Be as a function of 
neutron energy. 

very heavy systems the Coulomb force becomes an impor- 
tant effect and manifests itself as neutron excess (or proton 
deficiency) for the stability of a nucleus. We found that it 
is necessary to modify the global value of D (Eq. 13) for 
nuclei heavier than lead to account for this effect. We 
determine that the best value of D for, Z T > 82 is given 
by 

D = D g — £ t / (A t — Z t ). (16) 

For any regions not covered by these special effects the 
global value of D (Eq. (13)) gives good results. 

3. Results /conclusions 

Figs. 1-12 show the plot of the total absorption cross 
sections for the neutron-nucleus collisions. We observe 
that cross sections are never negative and are always 
positive quantities. For a particular system the energy at 
which the cross section drops down to zero is the lowest 
energy to which our model should be used. The data for 



Fig. 2. Total absorption cross sections for n + ^ 2 C as a function of 
neutron energy. 



Fig. 3. Total absorption cross sections for n + ^Al as a function of 
neutron energy. 


higher energies (379-1731 MeV) has been taken from Ref. 
[10], and for lower energies has been taken from the 
compilations of ENDF/B-VI, and Refs. [1 l,12].The agree- 
ment with experiment is excellent for all nuclei for the 



Fig. 4. Total absorption cross sections for n -F^jFe as a function of 
neutron energy. 



Fig. 5. Total absorption £ross sections for n 4- 28 Ni as a function of 
neutron energy. 
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Fig. 6. Total absorption cross sections for n + ^Cu as a function 
of neutron energy. 


entire energy range. We have also compared (Figs. 2, 3, 6, 
8, 10) present results with the two popular models: model 
developed by Wilson et al. [6] here at the NASA Langley 
Research Center (LaRC (Prior)), and that of Letaw et al. 



Fig. 7. Total absorption cross sections for n + 5£Zn as a function 
of neutron energy. 



Fig. 8. Total absorption cross sections for n + J? 7 Ag as a function 
of neutron energy. 



Fig. 9. Total absorption cross sections for n + ^Sn as a function 
of neutron energy. 

[7] at the Naval Research Laboratory (NRL). We see from 
the plots that both these models agree very well with 
experiment at higher energies, however, there are depar- 



Fig. 10. Total absorption cross sections for n + g 2 8 Pb as a function 
of neutron energy. 



Fig. 1 1. Total absorption cross sections for n + g^Bi as a function 
of neutron energy. 
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Fig. 12. Total absorption cross sections for n + as a function 
of neutron energy. 

tures at lower energies. In general LaRC gives better 
results than NRL model at lower energies. 

Clearly present model (LaRC) improves on earlier 
models and gives accurate results to much lower energies, 
and is a welcome improvement for the transport model 
data base. Neutrons are the dominant component of the 
radiation environment for the High Speed Civil Transport 
(HSCT) mission where the data base proposed here is of 
immense importance. 

We have successfully presented simple accurate and 
unified model for the total absorption cross sections for 
any system of charged and/or uncharged collision systems 
valid for the entire energy range, which are important 
ingredients for several data bases including transport. The 
model will also be useful where analytical forms of the 
cross sections are needed. 

In view of its utility the computer program calculating 
the comprehensive accurate absorption cross sections of 
paper I and present paper will be placed in the COSMIC 
program library for the public use. 


Note added in proof 

There is an error in Eq. (12) of paper I (Ref. [4]). A 
pair of parentheses is missing. The correct equation should 
read as follows: 

D -2.77 — 8.0 X 10“ 3 A t + 1 .8 X 10" 5 A^ 

- 0.8/( 1 4- exp((250 - £)/75)) . (12) 
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